A model for hand-over-hand motion of molecular motors 
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A simple flashing ratchet model in two dimensions is proposed to simulate the hand-over-hand 
motion of two head molecular motors like kinesin. Extensive Langevin simulations of the model 
are performed. Good qualitative agreement with the expected behavior is observed. We discuss 
different regimes of motion and efficiency depending of model parameters. 



I. INTRODUCTION 

The dynamics of molecular motors is an important 
topic in biophysics and nanotechnology. In the living 
and in the artificial nanoscale world fast non-diffusive di- 
rected transport or rotary motion constitute key ingredi- 
ents of any complex structure. Molecular motors are the 
"nanomachines" which perform these tasks [2, EL HI 0] • 
This definition involves a considerable amount of differ- 
ent molecules: Motor proteins, such as myosin and ki- 
nesin, RNA polymerases, topoisomerases, ... 

In this paper we focus on the problem of directed mo- 
tion over a substrate which is exemplified in the kinesin 
(E) @- Active transport in eukariotic cells is driven by 
complex proteins like kinesin which moves cargo inside 
cells away from the nucleus along microtubules trans- 
forming chemical fuel (ATP molecule) into mechanical 
work. Kinesin is a two head protein linked by a domain 
(neck) and a tail which attach a cargo or vesicle to be 
carried. The two heads perform a processive walk over 
the substrate (the microtubule). The way in which this 
process is performed attracts big interest in the research 
in molecular biology as well as in biological physics. In 
order to understand how kinesin works two properties 
that arise from the structure 0, 0] of the microtubules 
cannot be forgotten: they have a regular, periodic struc- 
ture and structural polarity - they are asymmetric with 
respect to their two ends, which determines the direction 
of kinesin motion 

In the last fifteen years, experimental molecular bi- 
ology has provide a lot of new results which allows to 
elucidate, at mesoscopic level, the main mechanisms for 
directed transport. These experimental evidences are 
mostly based in single molecules experiments [8, 9]. The 
interpretation of these results is not always easy and 
many times are not conclusive on the detailed way in 
which the motor walks. Two basic mechanisms have been 
proposed to explain the kinesin motion: "inchworm" and 
"hand-over-hand" motion (see figure[T]). In the first case, 
one head does not overcome the other one. In this case 
the period of the motion is one period of the microtubule 
structure (Iq in the figure) . In the hand-over-hand mech- 
anism one head overcomes the other. Now the period 
for each head is the double (21q). In both cases the 



Hand-over-hand Inchworm 
■A ft :» ™* i ft ft ;» m& 




FIG. I: (Color online) Schematic representation of the pos- 
sible mechanisms of motion for the motor of kinesin. In the 
hand-over-hand case each head moves a distance equals to 2lo 
whereas in the inchworm the period of motion is Oo- 



center of masses advances the same length. Although 
first single molecule experiments were compatible with 
both mechanisms more recent experiments have shown 
[Tol [ill . H2I [T3I ] in a very clever way that hand-over- hand 
motion may be more plausible. 

Two strategies can be devised in order to model the 
motion of molecular motors [3, HH ■ On one hand con- 
tinuous models based on mirror symmetry breaking po- 
tentials (ratchet potentials [HI]) or time symmetry broken 
driven forces [IT]]- On the other hand, discrete kinetics 
models which are based on the solution of master equa- 
tions associated to different states of the motor (see [14| 
and references therein). Using either approximation both 
mechanisms have been studied: inchworm [IH [l9j or 
hand-over-hand @, EH, El El . 

In this work we study a minimalist mechanical contin- 
uous model for hand-over-hand motion that, we believe, 
captures the main features of biological motors. The 
model also has into account the properties of the mi- 
crotubule substrate. The article is organised as follows: 
first we cast a two dimensional model which can mimic 
the motion of the motor. Within of reasonable values 
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of parameters we explore different regimes of motion. In 
the conclusions section we will discuss the validity of the 
results to model a molecular motor. 



II. 2-D MODEL 

In order to find a suitable model for the kinesin motor, 
its properties shall be studied carefully. Ref. @ sum- 
marizes all these features: kinesin is a two-head protein 
which moves along the microtubule with 8.3 nm steps, 
matching the repeat distance of the microtubule lattice; 
each step needs 1 ATP which is hydrolyzed and the move- 
ment stalls when a backward load of 7pN is applied. 
Experiments reported in [101 ] show, by marking one of 
the heads, that the motion follows the hand- over- hand 
mechanism, as a 16.6 nm step is observed for each head, 
thus forbidding the movement proposed in the inchworm 
mechanism. 

The description of the movement is rather simple: the 
two heads of the kinesin are attached to the microtubule 
[24| in two neighbor monomers until 1 ATP molecule is 
hydrolyzed by the head backwards. This energy frees 
the head, which moves to a new binding place ahead the 
other one. Two complementary mechanisms to under- 
stand how the particle released is able to find the next 
binding site have been proposed 0, 12a ]: (a) The neck 
linker mechanism assumes a conformational change in 
the neck between heads which moves the free head from 
one place to the next forwards, (b) The diffusional search 
relies on the assumption that the noise associated to the 
thermal bath that surrounds the particle makes the free 
particle move, and this movement is preferably forwards 
and forced by the particle ahead which is attached to the 
microtubule. 

Thermal fluctuations play a central role in the whole 
process. In the nanometer-length dimension and at room 
temperature, motion is governed by randomness induced 
by the environment (in this case the cytosol, made up 
mainly by water). At this scale damping and thermal 
noise are dominant and the dynamics can be studied by 
an overdamped Langevin equation: 
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dr 
~dt 



-W(r)+F(r, t)+£(t). 



(1) 



Here F stands for external forces and £ for thermal noise, 
being 

<&(*)&(*') >=2~/k B T6(t-t')6 jk (2) 
(£j and £fc are Cartesian components of the vector £). 

A. Energy potentials 

We will model the kinesin as two interacting particles 
moving in the plane under the effect of flashing ratchet 
substrate potentials (two particles moving in two dimen- 



The potential energy of the system is given by 
V(ri,r a ) = Vi(ri,t) + V 2 (r 2 ,t) + V 12 (r 1 - r 2 ) 



(3) 



The two heads of the kinesin are linked through a mod- 
ified version of the Finite Extensible Non-linear Elastic 
(FENE) interaction [26|: 
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with r = \tx — ral, K is the stiffness of the neck, Iq the 
equilibrium distance between heads and Rq determines a 
maximum allowed separation, Iq — Rq < r < Iq + Rq. 

With respect to the substrate potentials, in order to 
model the characteristics observed, two periodic flashing 
ratchet potentials lagged half a period in the x direction 
will be used. 



V j (r,t) = V j (r)f j (t). 



(5) 



j = 1,2 and in the x direction the potentials are periodic 
with period 21q and Vo is displaced Iq, the period of the 
microtubule lattice [27[, with respect to V\\ 



V 1 (r + 2l Sc)=V 1 (r)=V 2 (r + l x) 



(6) 



The mathematical description of the 2d potential as- 
sociated to particle 1 [see Fig. ([2])] is the following 



V 1 (x,y) = V lx (x) + V ly (y) 



with 



V lx (x) 



: Vq if < x < Xm 



iBtV if x M <x< 2l . 



(7) 



(8) 



xm controls the asymmetry of the potential and if xm — 
Iq the potential is symmetric. 

In order to confine the particles in the microtubule 
channel, with respect to the y direction we choose a sim- 
ple parabolic dependence. 
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We still have to define fjit). The idea is to reproduce 
a cyclic motion. Such a cycle has 4 steps, see Fig. ([31 
First (i = 0) both particles are confined close to the 
minima of their respective potentials and thus separated 
an averaged distance Iq (the natural length of the neck). 
After a given time t on some energy arrives at particle 1 
for instance which does not see its substrate potential for 
a time t a Q. During this time this particle suffer a thermal 
diffusion only subjected to the interaction with the other 
particle. When V\ is switched on again at t = t on + t ff , 
the particle slides down towards some minimum energy 
position. This step lasts another t on time and then at 
t = 2t on + i Q ff , Vz is switched off for a t a g time closing 
the cycle. The total period of this cycle is T = 2t on +2t g. 
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FIG. 2: (Color online) Surface plot of the 2d substrate po- 
tential Vi- Minima correspond to x — 2n (n = 0, ±1, ...) and 

y = o. 
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FIG. 3: (Color online) Time sequence for the flashing ratchet 
substrate potentials. Each potential acts on a different par- 
ticle. Note that this potential follows the sequence of attach- 
deattach shown in the hand-over-hand motion, Fig. [T] 



B. Normalization 

We will measure distance in units of Iq = 8.3 nm, 
the distance between monomers in the microtubule, see 
also 27]. Energy is measured in units of Vo, the max- 
imum value of the substrate potential. We choose 
V ~ E ATP - 20 k B T (at 300 K) ;28]. The natural 
unit of time will be r = IqJ/V) ~ 40 ns. Here, 7 is 
the damping coefficient used in the Langevin equation 
(7 = Qm]r = 4.7 • 10~ n ^, with 77 = 10~ 3 Pa the viscos- 
ity of the water and r = 25 A the size of the head) . 

We will use now ~ signs for normalized variables as 
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RESULTS 



We are going to present our results based in the nu- 
merical integration of the normalised system of equations 
for the two particles. The integration algorithm we use 
is a version of the Runge-Kutta algorithm for integration 
of stochastic differential equations (3o4s2c) [29|, |30| . 
With respect to the different constants and parameters, 
unless extra information is given, the default normalized 
parameters will be we T — 0.05 (300 K), i ff = t on = 20, 
K = 10, k y = 1, Ro = 0.4 and x M = 0.5. 



A. Dynamics of the system 



As we will see, thanks to the asymmetric character of the 
potential a directed motion is obtained. 

In order to compute the efficiency of the motion we 
define an efficiency parameter given by: 



(Asi) 
2Z 



x 100, 



(10) 



where (Axi) is the average advance of particle 1 (for 
instance) per cycle of the potential. Note that our def- 
inition of efficiency is basically the velocity of the mo- 
tion, in fact the mean velocity can be computed as 
Wmcan = x 2lo/{t on + t„f f) and is not related to the 
input of energy and the output of work. 

The important parameter here is t s, the time a par- 
ticle has for the diffusive motion. t on only requires to be 
long enough to allow for relaxation towards a minimum, 
which in over damped dynamics happens very fast. Thus, 
in our simulations we have played with different values 
of t g and set t on = t g. This value corresponds to a 
duty ratio r = t on /(t on + t ff) = 0.5 which guarantees the 
processivity of the motion [J, [ll| . 



Fig. [4] shows a typical example of the dynamics of the 
system at the parameter values listed above. There we 
can see that simulations reproduce the expected mech- 
anism, a hand-over-hand net advance of the molecule. 
Middle figure shows a detail of the top one. 

If both potentials are on, the particles do random 
motions around the minimum potential energy position. 
However, as one of the potentials is turned off, its linked 
particle starts to diffuse in 2d. The importance of the 
asymmetric mechanism is fully understood here. After 
t ff, when the potential is turned on again, most of the 
times the particle is sited to the right of the maximum of 
the asymmetric potential and then typically moves down 
to the nearest minimum position. As we have said, due 
to the asymmetry of the potential, this minimum more 
frequently corresponds to the one to the right of the orig- 
inal one. Clearly, the more asymmetric the potential is, 
the more likely the system moves forward. 

In the third graph of Fig. [4] we show the trajectories 
of the particles in phase space. The distance between 
heads moves around the rest distance lo- Motion in the 
x direction happens usually when one of the substrate 
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FIG. 4: x(t) for the two particles (top and middle) and tra- 
jectory in the x-y plane (bottom). The middle figure also 
shows the flashing dynamics of the substrate potentials (the 
base lines correspond to the on periods) . 



FIG. 5: x-axis projection, at different normalized times t, of 
the diffusion of a particle attached to the other, fixed at (1,0), 
when no substrate potential is being applied. Data obtained 
at f = 0.05 



potentials is off. Otherwise particles stay most of the 
time close to minimum energy position. 



B. Efficiency as a function of t s and T 

A hrst estimation of the time needed for a particle to 
reach the next minimum can be easily worked out by us- 
ing the 2D diffusion equation for the particle probability 
distribution p(ip, t), 



(12) 



Writing |12| in polar coordinates and assuming that 
the distance between heads r is constant the equation 
reads 



dp 

~dt 



D 



d 2 p 1 dp d 2 p 
dr 2 r dr dip 2 



D d 2 P 
r 2 dip 2 



(13) 



This equation can be solved (making Fourier transforma- 
tion for instance) with appropriated initial conditions 



to give the normalized p(ip, t) 
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(15) 



from this results the mean angle reached at time t is given 
by 



<^ 2 > 
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(16) 
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FIG. 6: Efficiency as a function of i ff for different normalized 
temperatures T. 



FIG. 7: Efficiency as a function of t s at different positions 
of the maximum xm- When xm = 1.0, the potential becomes 
symmetric and no rectified movement is observed. 



where we have used the Stokes-Einstein relation: 



D 
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(17) 



Let be y>M the angle where the maximum of the po- 
tential is placed which is determined by the x position 
of that maximum, xm- If the particle is at a i posi- 
tion less than xm when the potential turns on, it will 
return to its original position. However, if x > xm the 
particle will move forward. Assuming that the position 
of the maximum is placed at Xm = 0.5, and r — 1 we 
obtain ipM — 7r/3. Figure [5] shows the time evolution 
of the probability distribution (projected on the x-axis). 
It is clearly observed that as time goes the probability 
of finding that a particle crosses the maximum xm in- 
creases. The time in which the probability for crossing is 
1/2 is simply given by 



t = 0.674 



2k R T ■ 



(18) 



With the values given above the adimensional time (for 
temperature T = 0.05) is t ~ 16. For this t s time the 
efficiency of the motor is half of the maximum one, which 
is fixed by xm (see below). 

Finally we analyze the behavior of the efficiency with 
temperature, Fig. [6l For low temperatures we need long 
t a ff times to reach a reasonable efficiency as expected 
from equation |18| and we not get an asymptotic limit. 
For intermediate temperatures T = 0.03 — 0.05 the high- 
est efficiency is achieved. Moreover, in the limit of high 
temperatures compared to the 2d potential and long t a s, 
the efficiency starts to fall, as backwards movement is 
more likely to occur (the free particle can drag the con- 
fined one). 



C. Efficiency at different asymmetries 

Here we present results on the behavior of the sys- 
tem as the asymmetry of the potential changes, being 
xm the parameter that controls it (x m fixes the position 
of the maximum in the period 21q periodic potential, so 
%m = 1 corresponds to the symmetric case). At a given 
t ff time, efficiency depends importantly on this param- 
eter. The mechanism will be inefficient for a symmetric 
potential and the largest efficiency will be obtain for the 
more asymmetric one. 

Fig. [3 shows the numerical simulation of the efficiency 
as a function of f ff f° r different values of xm- As we 
reduce the asymmetry the efficiency tends to zero, as 
shown in the Xm = 1-0 line, which corresponds to a sym- 
metric potential. On the other hand the efficiency of 
the mechanism increases as we make the potential more 
asymmetric. In all the cases, when we increase t s the 
efficiency grows from zero and saturates at its maximum 
value for long enough values of this parameter. 



D. Dynamics under external loads 

In this section we want to explore the experimental 
results reported in Ref . || , where backward stepping was 
observed when using high backward loads. Then, it is 
worth studying how the system behaves under the effect 
of an external force. 

To model the effect of such a load is not trivial. We 
have to decide how the total load Q is divided into the 
two heads of the protein. It seems obvious that a head 
can make an opposite force to the applied only in case 
it is fixed to the microtubule. Therefore, the following 
mechanism is proposed: if there is only one head with 
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FIG. 8: Efficiency as a function of the external load applied, 
Q. Each line refers to a different t ff- 



FIG. 10: Efficiency as a function of the natural length of the 
neck, 1' , for different values of t f{. 
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FIG. 9: Efficiency as a function of the external load applied 
Q, for different normalized temperatures T. 



its potential switched on, it will bear the whole oppo- 
site load. On the other hand, if both heads have their 
potential on, everyone will bear a force Q/2. 

The expected behavior of the system is the following: 
as one potential turns off, its associated head starts dif- 
fusing. The other particle feels a force Q which dou- 
bles the previous Q/2. Therefore, if that force is strong 
enough, the particle starts climbing the potential slope. 
The asymmetric potential plays again an important role: 
if the external force is positive, the particle faces the 
sharpest slope of the potential, so a bigger force than 
in the negative case is needed. 

Fig.[5]shows the relationship between external load and 
t ff. The most important characteristic is the value of the 
load for which the system does not move, efficiency. For 



negative loads, as t g shortens, the particle needs greater 
forces to start to move backwards. On the other hand, 
when long times are employed, the mechanism seems to 
reach a limit around Q = —0.5. 

We have also studied the effect of the temperature in 
the mechanism. Results are shown in Fig. [5] where effi- 
ciency versus external load for a given value of i ff = 20 
is plotted at different temperatures. An almost linear re- 
lation between critical load and temperature is obtained 
in this range. 



E. Varying the natural length of the neck 

Up to now we have studied the case where the 
two space lengths of the system, the distance between 
monomers in the microtubule and the natural length of 
the neck, are equal (both are Iq). In this section we have 
extended our work to the study of the case when the 
natural distance between the heads is different from the 
spatial unit, fixed by the distance between monomers in 
the microtubule. Then, in our model, Iq need to be re- 
placed by 1' in Eq. g]). 

Fig. [10] is clear enough to provide strong evidence 
about the striking behavior observed as the natural 
length of the neck tends to 0: 100% efficiency is achieved. 
This almost deterministic mechanism can be understood 
with the help of Fig. [TT| and presents three steps: (a) 
We start with one particle sited in a minimum of the 
potential and the other one ahead (it feels a small force 
since the potential slope there is also small), (b) As the 
first potential disappears, the second particle moves to 
its minimum, dragging the other one. (c) Now the po- 
tential turns on, thus making the first particle to move 
ahead and we recover a situation equivalent to step (a). 

Fig. [TZl shows results at different temperatures. First 
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FIG. 11: (Color online) Schematic explanation of the almost 
deterministic motion observed when 1' close to zero. 
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FIG. 13: Efficiency versus 1' for different values of K. 

The study of the efficiency as a function of 1' at different 
values of K is shown in Fig. [T3j For l' Q < 0.5, the stiffness 
of the neck determines whether the particles prefer to be 
in their minimums no matter how far they are, or in an 
intermediate position, as plotted in Fig. 1111 a). 
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FIG. 12: Efficiency versus 1' for different temperatures, T. 



of all, the deterministic T = limit must be carefully 
explained. In this limit, there are only two possible 
values for the efficiency: 0%, associated to the range 
% E (0.5,1.5), and 100% for \ £ [0.0, 0.5) U (1.5,2.0]. 
These two regions can be fully explained using the mech- 
anism described above. In the 0% case, switching off one 
of the potentials will make the other particle move to a 
minimum, but it will not be the minimum ahead which 
would not produce a net movement forwards. 

There is just one parameter left to be discussed, which 
is the stiffness of the linker between heads of the motor. 



IV. CONCLUDING REMARKS 

We have studied a simple mechanical model for hand- 
over-hand motion in two dimensions. This model has into 
account some important characteristics of two heads bio- 
logical motor as kinesin. These characteristics are incor- 
porated in the model in a simple but realistic way. The 
hand-over-hand mechanism requires a two-dimensional 
space. Unidirectionality is given by the ratchet potential 
in the advance direction. The balance between on and 
off times controls the efficiency and processivity of the 
motion. With all these ingredients we have been able to 
simulate the most remarkable features of kinesin motion 
within reasonable values of he parameters. Specifically, 
we have clearly observed a stochastic directed motion in 
which particles alternates each other (hand- over- hand). 
Moreover, a strong dependence of the stall force with 
off time and temperature has been found. Temperature 
makes a decrease of stall force with respect to one ex- 
pected from energetic calculations. This decrease in the 
motor efficiency agrees with experimental observations 

am. 

Several improvements to the model can be considered 
in future work. A link between t Q g and ATP concentra- 
tion could be established. This would imply a random 
flashing force instead of the periodic one used here. An- 
other interesting extension of the model could allow the 
motor to change the lane in y axis. This could be easily 
implemented by using a periodic potential in the trans- 
verse direction. 
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Finally, we have to stress that the characterization 
of the behavior and properties of those motors and the 
mechanisms behind them is an initial step toward the 
construction of synthetic nanoscale motors. This is a 
very active field in the nanoscicncc world. There has 
been some successful achievements in this field that in- 
clude triptycene motors [32| . helicene motors [33TJ and a 
nanotube nanomotor (34j . In this article, we have shown 
the conditions for which a nanowalker can work. 
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